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ABSTRACT 

We present the results of a survey for Extremely Red Objects (EROs) undertaken 
in the fields of ten massive galaxy cluster lenses at z ~ 0.2, combining sensitive, 
high-resolution Hubble Space Telescope imaging with deep, half-arcsecond X-band 
imaging from UKIRT. We detect 60 EROs with (R - K) > 5.3, of which 26 have 
(R — K) > 6.0 in a total image plane survey area of 49arcmin 2 down to K = 20.6, 
including one multiply- imaged ERO. We use detailed models of the cluster lenses 
to quantify the lens amplification and thus correct the observed number counts and 
survey area for the effects of gravitational lensing. After making these corrections, 
we estimate surface densities at K < 21.6 of (2.5 ± 0.4) and (1.2 ± 0.3) arcmin~ 2 
for EROs with (i? — K) > 5.3 and 6.0 respectively. These ERO number counts agree 
with previous shallower surveys at K <, 19 and flatten significantly at magnitudes 
fainter than K ~ 19-20. This flattening may be due to a transition from an ERO 
population dominated by evolved galaxies at z ~ 1-2 (K <> 19.5) to one dominated 
by dusty starburst galaxies at z > 1 (K <; 19.5). We also compare our results with 
various model predictions, including a model that attempts to explain EROs in terms 
of a single population of elliptical galaxies formed at high redshift. We find that a 
formation epoch of Zf ~ 2.5 for this population matches the observed surface density 
of (R- K) > 5.3 EROs quite well, and the (R - K) > 6.0 sample less well. More 
sophisticated models, including semi-analytic prescriptions, under-predict the ERO 
surface density by approximately an order of magnitude, suggesting that these models 
produce insufficient stars and/or dust at high redshift. This deficit of EROs appears 
to be a general problem with models that reproduce the median redshift from K- 
selected redshift surveys. One possible explanation is that the current K- selected 
redshift distribution may be incomplete beyond z ~ 1. 
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1 INTRODUCTION 

Extremely Red Objects (EROs) are defined by their very 
red optical/near-infrared colours, e.g. (R — K) > 5.3 (e.g. 

* E-mail: graham.smith@durham.ac.uk 



Moriondo et al. 2000) and typically have Tf-band magni- 
tudes of K <; 17.5. Since their discovery (Elston et al. 1988, 
1989, 1991; McCarthy et al. 1992; Eisenhardt & Dickinson 
1992; Graham et al. 1994; Hu & Ridgway 1994; Dey et al. 
1995), they have intrigued cosmologists, but until recently 
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their extreme properties have escaped detailed, systematic 
investigation. 

The variety of colour criteria used to define EROs e.g. 
(R - K) > 5, > 5.3, > 6.0 or (I - K) > 4, > 5, together 
with different survey depths and apparently strong cluster- 
ing (Daddi et al. 2000a) have led to a wide dispersion in the 
claimed surface densities of EROs (e.g. Hu & Ridgway f994; 
Cowie et al. 1994; Aragon-Salamanca et al. 1994; Dey et al. 
1995; Moustakas et al. 1997; Yamada et al. 1997; Beckwith et 
al. 1998; Eisenhardt et al. 1998; Barger et al. 1999; Thomp- 
son et al. 1999; Daddi et al. 2000a; McCracken et al. 2000). 
The larger and newer surveys suggest that to K % 19 the sur- 
face density of (R — K) > 6.0 galaxies is 0.17±0.07 arcmin -2 
(Daddi et al. 2000a), increasing to 0.28 ± 0.04 arcmin ~ 2 at 
the faintest limit probed to date K 20 (Thompson et al. 
1999). 

The optical faintness of EROs (R ^ 24) and the lack 
of near-infrared multi-object spectrographs have combined 
to severely limit the number of spectroscopic identifications. 
However, the available spectra indicate that EROs comprise 
two broad classes of system: galaxies with evolved stellar 
populations at z ^ 1 (e.g. Dunlop et al. 1996; Soifer et 
al. 1999; Liu et al. 2000) and dust-reddened starbursts at 
similar and higher redshifts (Graham et al. 1996; Dey et 
al. 1999). A few starburst EROs have also been detected 
in radio and sub-millimetre (sub-mm) wavelength surveys 
of the distant universe (Spinrad et al. 1997; Smail et al. 
1999a, 2000; Dey et al. 1999; Gear et al. 2000; Barger et al. 
2001). However, morphological studies (e.g. Treu et al. 1999; 
Moriondo et al. 2000) suggest that at K ~ 19 the bulk of 
EROs are evolved elliptical galaxies. 

These two classes of EROs may represent different 
phases in the formation and evolution of a single family 
of high redshift galaxies: massive ellipticals. However, rival 
galaxy formation theories (pure luminosity evolution and hi- 
erarchical clustering) predict very different formation epochs 
for such galaxies. The former predicts formation at high red- 
shift in what is traditionally known as the "monolithic col- 
lapse" scenario (e.g. Eggen et al. 1962; Larson et al. 1975; 
Tinsley & Gunn 1976) followed by subsequent passive evo- 
lution. The latter predicts formation through the merging of 
disk galaxies at z ^ 1 (e.g. White & Frenk 1991; Kauffmann 
et al. 1993, 1996; Baugh et al. 1996; Cole et al. 1996, 2000). 
The properties of EROs, specifically, their redshift distribu- 
tion and spectral characteristics, are therefore a crucial test 
of galaxy formation theories and will improve our ability to 
make complete measurements of the star formation history 
of the universe (e.g. Madau et al. 1996; Smail et al. 2001). 

Many groups are currently investigating different as- 
pects of the ERO population using a variety of observational 
approaches, for example: wide-field imaging surveys (Daddi 
et al. 2000a; Thompson et al. 1999), morphological studies 
(Treu et al. 1999; Moriondo et al. 2000), searches for over- 
densities around radio-loud AGN (Willott et al. 2000; Chap- 
man et al. 2000; Cimatti et al. 2000), deep Tf-band number 
counts (McCracken et al. 2000; Barger et al. 1999), searches 
for sub-mm counterparts (Smail et al. 1999a), photometric 
redshift studies (Fontana et al. 1999), studies of faint field 
spheroidals (Menanteau et al. 1999) and high redshift hosts 
of optically faint X-ray sources (Cowie et al. 2001). 

This paper describes a survey whose aim is to produce 
a sample of EROs which is optimised for follow-up, with the 



specific goal of improving the opportunities for measuring 
redshifts and spectral properties. To achieve this we take 
advantage of the natural magnification of massive clusters of 
galaxies to amplify background source populations, in this 
case EROs, and so boost their apparent brightness. This 
approach has been successfully applied to the investigation 
of other faint sources with steep number counts, such as 
mid-infrared and sub-mm galaxies (e.g. Altieri et al. 1999; 
Smail, Ivison & Blain 1997) and may similarly aid progress 
in understanding the properties of EROs. 

We use a unique sample of X-ray luminous galaxy clus- 
ters at moderate redshifts which have been uniformly imaged 
with the Hubble Space Telescope (HST) through the F702W 
filter to provide deep, high resolution imaging of the clus- 
ter cores - including many gravitationally lensed arcs. We 
have supplemented these deep optical images with high qual- 
ity ground-based near-infrared imaging with UKIRT and we 
use this combined dataset to construct a sample of gravita- 
tionally amplified EROs. 

In the next section we describe the optical and near- 
infrared imaging and its reduction. §3 describes the analysis 
of these images and the properties of the ERO sample con- 
structed from them. Finally in §4 we discuss our results and 
list the main conclusions of this work. Throughout this pa- 
per we assume Ha = 50kms _1 Mpc _1 , £7o=0.3 and A = 0.7. 



2 OBSERVATIONS 

2.1 HST/ WFPC2 Optical Imaging 

The HST imaging of the clusters used in this paper was 
obtained as part of a survey of gravitational lensing by X-ray 
luminous clusters (Smith et al. 2001a, 2001c). The complete 
sample consists of twelve X-ray luminous clusters (Lx > 
8 x 10 44 ergs -1 , 0.1-2.4 keV) in a narrow redshift slice at 
z ~ 0.2, selected from the XBACs sample (Ebeling et al. 
1996). For the purposes of this paper we concentrate on the 
ten clusters listed in Table 1 for which homogeneous F702W 
imaging is available (the remaining two clusters which are 
excluded from our analysis are A 1689 and A 2390 for which 
only archival F814W observations are available). 

Of these ten clusters, A 2218 was observed in Cycle 4, 
A 2219 in Cycle 6 and the remaining eight clusters in Cy- 
cle 8. All clusters except A 2219 were observed for three, 
single orbit exposures with the WFPC2 camera through 
the F702W filter; A 2219 was observed for six orbits using 
the same instrument/filter combination. The total exposure 
times for the fields, T_r, are listed in Table 1. Each expo- 
sure was shifted relative to the others by ten WFC pixels 
(~ 1.0") providing a partial overlap of the WFC chips. Af- 
ter pipeline processing, standard iraf/stsdas routines were 
used to align and combine the frames to remove both cos- 
mic rays and hot pixels. Corrections for under-sampling of 
the point spread function (PSF) and geometric distortion of 
the optics were made using the dither package within iraf 
(Fruchter & Hook 1997). More details are given in Smith et 
al. (2001a, 2001c). 

We adopt the photometric system from Holtzman et 
al. (1995) and, in the following, have chosen to convert the 
-R702 photometry to Cousins 72-band to aid in comparing 
our results with previous surveys. To achieve this we adopt 
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Table 1. Details of F702W and K-band Observations Used to Construct the ERO Sample 
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a (V — R) colour for an Sbc galaxy at z = 1-2 of (V — R) ~ 
1.1, based on spectral templates derived from local galaxies. 
This translates into a correction of R — R702 ~ 0.4 with a 
systematic uncertainty of g 0.1 magnitudes arising from the 
likely presence of both more and less evolved galaxies than 
the adopted Sbc spectral type. The equivalent correction 
for galaxies with the (V — R) colours of cluster ellipticals at 
2 ~ 0.2, (V-R) ~ 0.9, is R-R 70 2 ~ 0.3, with an estimated 
systematic uncertainty of g 0.05. Both of these corrections 
are sufficiently accurate to define a sample of EROs in these 
fields. 

The final WFPC2 frames have an effective resolution 
of 0.15" and a typical 3-ct detection limit within our 2"- 
diameter photometry aperture of R ~ 26.6. 

2.2 Near-infrared Imaging 

The near-infrared (if-band) imaging essential for identify- 
ing EROs was obtained for nine of the ten clusters with 
the UFTI imager on the 3.8-m United Kingdom Infrared 
Telescope (UKIRT)*. UFTI incorporates a 1024 2 Hawaii-1 
detector providing a 92" field of view with 0.0908" pixel -1 
sampling - necessary to sample the best seeing provided by 
UKIRT, £ 0.2-0.3". The A"-band imaging of these fields 
was obtained in two observing runs in 2000 April 4-7 and 
September 26-29. Conditions during both observing runs 
were good, with reasonable seeing and transparency. Each 
WFPC2 field was covered in three pointings (one per WFC 
chip: WF2, WF3 and WF4), although the different roll- 
angles for the various WFPC2 observations means that a 
particular near-infrared exposure may overlap more than 
one WFC chip. The individual pointings usually consist of 
54 x 120-s exposures grouped into six sets of nine exposures, 
or occasionally 72 x 90-s exposures in eight sets of nine, each 
exposure being offset on a 3 x 3 grid with 10" spacing, and 
the origin of each set was randomly moved by ~ 3". We list 

* UKIRT is operated by the Joint Astronomy Centre on behalf 
of the Particle Physics and Astronomy Research Council of the 
United Kingdom. 



the total exposure times for the individual pointings (Tjf) 
and the seeing for each in Table 1. The median seeing for 
these 2-hr integrations is 0.5". 

The data were reduced using the dedicated UFTI sum- 
mit data pipeline (Currie et al. 1999) including dark sub- 
traction, fiat fielding using a local median sky, resampling 
onto an astrometric grid and combining to remove defects 
and cosmic ray events. A small amount of the data was re- 
reduced in Durham following a similar scheme and using 
standard iraf tasks. Calibration was obtained from UKIRT 
Faint Standards (Casali & Hawarden 1992) interspersed 
throughout the science observations. We estimate that our 
absolute calibration is good to <; 0.05 mag. A reddening cor- 
rection of E(R — K) — 2.31E(B — V), was also applied to 
the (R - K) colours following Schlegel et al. (1998). The 
final UFTI Tf-band images typically reach a median 80% 
completeness of K ~ 20.6 (Table 1 & §2.4). For comparison 
a passively evolving L* elliptical will have K ~ 18-19 at 
2 = 1. We show an example of one of our fields in Fig. 1. 

The tenth cluster, A 2218, was observed in commission- 
ing time with the new INGRID near-infrared imager (Pack- 
ham et al. 2001) on the 4.2-m William Herschel Telescope 
(WHT)t. INGRID comprises a 1024 2 Hawaii-1 detector giv- 
ing a 248" field of view with 0.242" pixel -1 sampling. These 
data comprise a 8.3 ks exposure of the whole WFPC2 field 
in A 2218 (Kneib et al. 1996) in the K s band. The seeing 
on the final stacked frame is 0.75" and the frame reaches 
an 80% completeness limit of A" = 21.5 for point sources. 
The reduction, calibration and analysis of these data are dis- 
cussed in more detail in Smail et al. (2001). For uniformity 
in our analysis we restrict ourselves to a limit of K < 20.6 
in A 2218. 

In addition to the Jf-band imaging described above, 
a sub-set of the cluster fields were also observed in the J- 
band. We observed A 773, A 963, A 1763 and A 1835 with 

t Based on observations made with the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto dc Astrofisica dc Canarias 
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UKIRT/UFTI on the nights of 2001 February 16 and April 
4 in non-photometric conditions and ~ 0.6" seeing, and A 68 
A 267 and A 2219 with WHT/INGRID on the nights of 2001 
June 27 and 2001 May 6, in ~ 0.8" seeing and photometric 
conditions (Table 2) . All of these observations were reduced 
in a similar manner to the if-band images. The UFTI obser- 
vations reach a 5-cr limit of J ~ 21, whilst INGRID frames 
reach 5-cr depths of J ~ 20-21. Finally, a J-band image of 
A 2218 is also available from the INGRID commissioning 
observations of that cluster, however, as we find no ERO 
sources brighter than K = 20.6 within that field, we do not 
discuss these data further. 



2.3 Ground-based Optical Imaging 

To study the spectral energy distributions of EROs between 
the R— and K- bands, we also make use of archival /-band 
imaging of the clusters in our sample. These observations 
come from a number of ground-based 4-m class telescopes 
and comprise 0.5-3.6 ks integration in 0.7-1.4" seeing (Ta- 
ble 2), reaching I ;> 23. We use panoramic imaging of A 68, 
A 209, A 267 and A 383 from the CFH12k camera on the 3.6- 
m Canada- France-Hawaii Telescope, $ taken on 1999 Novem- 
ber 10-14. More details of their reduction and analysis can 
be found in Smith et al. (2001a, 2001c) and Czoske et al. 
(2001). For A 773 and A 963 we use /-band images taken 
on 1994 December 8-9 with PFCCD of the 4.2-m William 
Herschel Telescope^ . Finally, we exploit imaging of A 1763, 
A 1835, A 2218 and A 2219 taken with the COSMIC imager 
spectrograph on the Hale 5-m^, more details of these obser- 
vations and details of their reduction can be found in Smail 
et al. (1998) and Ziegler et al. (2001). 



3 ANALYSIS AND RESULTS 

3.1 Source Detection and Photometry 

The near-infrared images of each cluster were aligned and 
mosaiced together (if necessary) to create one fi'-band frame 
per cluster that covers the relevant WFPC2 field of view. 
The WFPC2 frames were rotated and aligned to the .fir- 
band frames with an rms tolerance of <; 0.01" and the whole 
field astrometrically calibrated using the APM catalogue to 
an absolute accuracy of 0.4". 

To produce a catalogue of EROs in these fields we first 
analysed the fi'-band frames using the SExtractor pack- 
age (Bertin & Arnouts 1996). All objects with isophotal 
areas in excess of 10 pixels (0.082 arcsec 2 ) at the fix = 
23 mag arcsec" 2 isophote (1.5<r pixel -1 ) and lying within the 

t Based on observations made with the Canada France Hawaii 
Telescope, operated on the island of Hawaii by the Canada France 
Hawaii Telescope Corporation 

§ Based on observations made with the William Herschel Tele- 
scope operated on the island of La Palma by the Isaac Newton 
Group in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias 

^ This work is based on observations obtained at Palomar Obser- 
vatory, which is owned and operated by the California Institute 
of Technology 



Table 2. Summary of /- and J-band Observations 
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(a) The /—band observations cover the whole WFPC2 field 
of view. For the J-band: WFPC2 denotes WHT/INGRID 
observations which cover the whole field; WF2/3/4 indicate 
UKIRT/UFTI imaging of a particular WFC chip. 

WFPC2 field of view were selected. Across the ten clus- 
ters, this survey covers a total area of 49 arcmin 2 , excluding 
the PC chips, and these catalogues contain a total of 2,382 
sources. We adopt the mag_best magnitude computed by 
SExtractor as the total fi'-band magnitude of each source. 

We perform extensive (~ 10 4 realisations) Monte Carlo 
simulations to measure the completeness limits of these 
catalogues by suitably scaling and re-inserting a moder- 
ately bright (K ~ 19), compact ERO source into the sci- 
ence frames. The resulting 80% completeness limits (roughly 
equivalent to 5-cr detection limits) are presented in Table 1. 
We also re-perform these simulations using a more diffuse 
galaxy of similar magnitude, and find that the typical 5-cr 
limiting depth was <; 0.3 magnitudes brighter than that for 
the compact source. This difference is not large enough to 
have a significant impact on the results presented in this 
paper. 

We measured the (R — K) colour of all these sources 
using a 2" diameter aperture on seeing matched R— and 
fi'-band frames. As the image smoothing required to match 
the seeing reduces the pixel to pixel variation in the sky 
backgrounds, the noise estimates from the seeing-matched 
frames are compromised. We therefore exploited the overlap 
regions between the three UFTI pointings used to cover each 
HST field to compare independent photometry of sources 
as a function of magnitude, and so derive conservative es- 
timates of the photometric uncertainties as a function of 
magnitude. We applied the same principle in the Zi-band, 
taking advantage of the six 2.4 ks exposures of A 2219 (§2.1), 
to make independent measurements on two frames, each de- 
rived from an independent sub-set of three exposures. This 
provides a reliable and conservative estimate of the photo- 
metric errors of sources in the /?-band. We also estimate the 
5-cr depth of the Zi-band imaging to be R = 26.0. However, 
for the purposes of our final photometry we adopt a 3-cr 
/?-band detection limit of R = 26.6 to discriminate between 
(R - K) > 5.3 and (R-K)> 6.0 (§2.4) sources lying close 
to both detection limits. 

We show the (R-K)-K colour-magnitude diagrams for 
the ten clusters in Fig. 2. The strong sequences of red galax- 
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Figure 1. An example cluster field from our HST / UKIRT survey, showing a true colour RIK image of the core of A 963. The circles 
indicate the thirteen EROs with (R — K) > 5.3 detected in this field. The two EROs north-west of the cD appear to be elliptical galaxies 
with similar redshifts on the basis of their optical/near-infrared colours and morphology (§3). The field is centred on the cluster cD and 
has North top and East left. Note the very blue giant arc visible south of the cD. 
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ies identifiable in all the panels correspond to the evolved 
early-type cluster members. These provide a useful check of 
the calibration of the (R — K) colours for sources in these 
fields. We estimate that the rms scatter in the colours of 
the cluster early-type sequences at a fixed luminosity as a 
function of redshift is ~ 0.08 mags, confirming that any field- 
to-field offset between our photometry of galaxies detected 
in each cluster field is < 0.08 magnitudes. 

Finally, we measure the (7 — K ) and (J — K) (where 
J-band imaging is available) colours of the ERO sample de- 
fined in §3.2 using a 2" diameter aperture on seeing-matched 
frames, adopting a 3-cr detection limit when no 7- or J- 
band counterpart is detected. 



3.2 ERO Selection 

We adopt the definition of an ERO as a galaxy with (R — 
K) > 5.3 as used by Daddi et al. (2000a), and also employ a 
more stringent definition of (R — K ) > 6 (Thompson et al. 
1999). These sample boundaries are shown for the individual 
cluster fields in Fig. 2, we list the numbers of galaxies in each 
class in the various clusters in Table 1 and give a catalogue 
of the sources in Table 3. 

We show in Fig. 3 the composite (R-K)-K colour mag- 
nitude diagram for all ten clusters, along with the sample 
boundaries. In total we find 60 sources with (R—K) > 5.3 in 
the ten fields (an image plane survey area of 49 arcmin 2 ) and 
26 with (R — K) > 6 down to our median 80% completeness 
limit of K = 20.6. 

In addition to the Monte Carlo simulations described 
in §3.1, we performed two further checks to verify the com- 
pleteness of our ERO sample. First, we visually checked 
true-colour optical/near- infrared images of the clusters (e.g. 
Fig. 1) to guard against losing ERO candidates due to 
contamination of photometric apertures by nearby sources 
with less extreme colours. In two cases we re-calculated 
the colours of objects after masking the light from nearby 
galaxies to ensure a more reliable measurement. Second, we 
searched for EROs hidden under the halos of bright cluster 
ellipticals by subtracting a median-smoothed version of the 
each Tf-band frame from the original frame to produce a 
"difference" frame. This search revealed no further EROs. 

These "difference" frames contain residual flux from the 
central regions of the bright cluster ellipticals, with the re- 
sult that a fraction of each WFPC2 field of view remains 
obscured. We estimate this fraction, first thresholding the 
"difference" frames at a level where the fainter EROs in 
our sample would not be detected. We then ray trace these 
thresholded frames back to the source plane (using the lens 
models described in §3.3). As the lens amplification is high- 
est at the cluster centre and the centre of other bright cluster 
ellipticals, the transformation to the source plane results in 
the obscured fraction of each field of view being typically 
£5%. 

We estimate contamination of the sample by red, low 
mass stars by visually comparing the Tf-band morpholo- 
gies of the 19 EROs in our sample with a FWHM < 
0.8" to the images of morphologically-classified stars se- 
lected from the HST frames. Only one of these sources 
(ERO J024805— 0330.2) has a star-like morphology; we 
therefore flag it as a possible star or AGN, but choose not 



to remove it from our analysis. This very low level of stellar 
contamination (< 1.5%) is consistent with that determined 
by Daddi et al. (2000a) and Thompson et al. (1999). 

A further potential bias in colour-selected surveys in 
gravitational lens fields is differential amplification across a 
galaxy image, causing separate regions of the galaxy with 
possibly different underlying colours to suffer different de- 
grees of amplification. As the lens amplification only varies 
strongly with position close to critical lines, this bias is only 
a concern for multiply-imaged galaxies. Only one of our sam- 
ple of 60 EROs is multiply imaged (§3.7). We therefore con- 
clude that this bias has a negligible effect on our sample. 

3.3 Surface Density of EROs 

In a "blank" field survey, the surface density of sources 
brighter than a given limiting magnitude can be calculated 
by dividing the number of detected sources by the surveyed 
area. This simple calculation is more complicated in the field 
of a gravitational lens, due to the amplification of the source 
flux and the accompanying distortion of the background sky 
(source plane). Consequently, some regions of the source 
plane are observed to greater depths than others, even if the 
image plane is observed to a uniform depth. These effects 
can be quantified by constructing a detailed model of the 
gravitational lens. Such models are then used to compute 
the lens amplification as a function of image plane position. 
This knowledge of the lens amplification allows the surface 
density of sources to be reliably estimated as described in 
more detail below. First, we briefly discuss the lens models 
used to quantify the lens amplification in our cluster sample. 

Smith et al. (2001c; hereafter S01) have constructed a 
detailed model of each cluster in our sample using the para- 
metric lens inversion method of Kneib et al. (1996). S01 
discuss the robustness of these lens models in considerable 
detail. In this paper, we summarise the redshift information 
used to constrain the models and thus provide an overview 
of the reliability of these models. The most robust mod- 
els are those constrained by the spectroscopic redshift of a 
multiply-imaged background galaxy. Four of the ten clus- 
ters fall into this category. Half of the remaining six clus- 
ters contain faint multiply-imaged sources, whilst only three 
contain no obvious multiple images in our F702W frames. 
S01 therefore used both photometric redshift techniques and 
uniformly selected samples of weakly sheared background 
galaxies to constrain the lens models for these six clusters. 
Consequently, the amplifying power of the four spectroscop- 
ically constrained models is known to ~2-5%, whilst the 
power of the other six is typically known to an accuracy of 
~ 10%. We therefore use SOl's lens models to remove the 
effects of lens amplification from both the source counts and 
the surveyed area using the method described below (Blain 
et al. 1999). 

The amplification suffered by each background galaxy 
depends on its redshift (z s ) and the redshift of the interven- 
ing lens (z L ). This redshift dependence is weak if z s ^> z L , 
and as the cluster lenses all lie at z L ~ 0.2 and the back- 
ground galaxies are all expected to lie at z s ~ 1-2, this 
regime applies. We therefore adopt a single source plane of 
z s =1.5 and use each model to compute a map of lens ampli- 
fication as a function of image plane position for this value 
of z s in each of the ten cluster fields. Adoption of a single 
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Figure 2. (il — K)—K colour- magnitude diagrams for the ten clusters in our survey, showing the order of magnitude field to field 
variation in ERO number counts. The _R"-band magnitudes in this figure are observed magnitudes (i.e. not corrected for gravitational 
amplification) and arc plotted down to the 80% completeness limits listed in Table 1 (median 80% completeness limit is K = 20.6). The 
two dashed lines on each panel indicate the (R — K) > 5.3 and 6.0 ERO selection criteria respectively. Points marked by an upward 
pointing arrow arc 3-cr limits in (R — K). 
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Figure 3. a) Composite (R — K)—K colour-magnitude diagram for the ten cluster fields in our survey, b) Zoom into (R — K) > 4 region 
of the composite (R — K)-K colour-magnitude diagram shown in (a). The two dashed lines in each panel indicate the (R — K) > 5.3 and 
6.0 ERO selection criteria. Galaxies marked by an upward pointing arrow are 3-cr limits in (R — K). We have co-added these objects and 
measure the typical i?-band magnitude (in a 2" diameter aperture) to be R = 27.0 ±0.5. This suggests that the galaxies not detected in 
the i?-band are a continuation of the (R — K) > 6.0 population, rather than a distinct population of galaxies with much more extreme 
colours. Note, the if -band magnitudes plotted in this figure are observed magnitudes (i.e. not corrected for gravitational amplification). 



value of z s introduces an uncertainty of ~ 10-20% into the 
final surface density values, which is comparable with the 
Poisson noise in the raw number counts. 

We first use the amplification maps to de-amplify the 
image plane flux of each ERO in our sample and hence 
obtain their source plane K -band magnitudes. As gravita- 
tional lensing is achromatic, no correction is required to the 
(R — K) colour. The number of EROs that are brighter than 
a source plane limiting magnitude ifn m , iV raw (< ifn m ), can 
then be found by simply counting the number of sources 
brighter than Ku m in the source plane after correcting for 
lensing. A simple Poisson uncertainty is attached to this 
value. 

We then calculate the area of the background sky within 
which each ERO is detectable in the relevant cluster field in 
the following manner. An ERO with a source plane mag- 
nitude of -Ksourcc will appear in the image plane of the 
relevant cluster with a magnitude brighter than an image 
plane detection limit of K dct if it is magnified by a fac- 
tor greater than ^ min = icr ' 4 ^^-^ "™'. The area in 
the source plane within which such a galaxy would be de- 
tected in that cluster is thus J 4>(/i m i n ), where Ay is the 
area of the z SOU rco = 1.5 source plane behind each cluster 
that lies within the WFPC2 field of view and is magnified 
by a factor greater than /i m in (Fig. 4) . Finally, we divide the 
number of galaxies iV raw by the sum of the areas (/^min) 
for all ten clusters, to obtain the cumulative surface density 

N(< Kum) =■ iVraw(< K hm ) /EAy (fl min ) . 

We present in Fig. 5 the cumulative surface density, 
N(< K), of (R-K) > 5.3 and (R-K) > 6.0 EROs detected 
in our survey after correcting the source fluxes and source 
plane surface areas for lens amplification. We also show 



the results of recent shallower, wide-field surveys (Daddi 
et al. 2000a; Thompson et al. 1999). The three datasets 
agree in the region of overlap (K ~ 17.5-19.5). We esti- 
mate the cumulative surface density of EROs at K <; 21.6 
with (R-K) > 5.3 to be (2.5 ± 0.4) arcmin~ 2 and with 
(R-K) > 6.0 to be (1.2 ± 0.3) arcmin" 2 . We note that 
the (R — K) > 6.0 EROs appear to comprise a constant 
fraction of the the overall (R — K) > 5.3 population at all 
magnitudes: ~ 0.40 ± 0.08. 

The slope of Daddi et al.'s (2000a) cumulative number 
counts for EROs with (R — K) > 5.3 is a = 1.05 ± 0.05 at 
K < 19.5, where N(< K) = U) aK . However, fainter than 
K ~ 19.5, we estimate the slope of our (R — K) > 5.3 
cumulative number counts to be a = 0.30 ±0.01, suggesting 
a break in the surface density of (R — K) > 5.3 EROs at 
K ~ 19-20. Number counts of the (R-K) > 6.0 population 
reveal a similar break, with slopes of a — 1.43 ± 0.25 and 
a = 0.37 ± 0.02 for brighter and fainter EROs respectively. 



3.4 Colour Distributions 

We present the (R — K)-(I — K) colour-colour diagram of 
our ERO sample in Fig. 6a. We also plot in Fig. 6a the ex- 
pected colours of passively evolving ellipticals as a function 
of redshift (for a formation redshift of Zf = 2.5 - see also 
§4.1). We begin by comparing the (R-K) and (I-K) ERO 
selection criteria. Of the (R - K) > 5.3 EROs, g 10 (15%) 
have (I — K) < 4 and would therefore not be classified as 
EROs in an (I — K) selected sample. Most of these galax- 
ies are also blue in (R — K), and given their proximity to 
the 2 <; 1 portion of the evolutionary track, we suggest that 
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Figure 4. Cumulative area of the source plane in our survey, 
A>, at z = 1.5, that experiences magnification greater than fi. 
The similarity of the individual clusters in this figure reflect the 
uniform selection criteria for our cluster sample (§2.1). 

they are probably early-type galaxies whose 4000-A break 
lies between the R- and /-bands. 

We also see that a substantial fraction of the (R — K) > 
6.0 galaxies are redder than the passive evolutionary track. 
Whilst we treat this comparison with caution, due to uncer- 
tainties in both photometric measurements and calibration 
onto the models, this suggests that the (R — K) > 6.0 popu- 
lation may contain a sizeable fraction of starburst galaxies. 
We also measure the typical colour of those EROs for which 
we can only place a limit on their (R — K) colour by co- 
adding the F702W-band images of these EROs. Using the 
same photometric approach as in §3.1, we obtain a typical 
7?-band magnitude for these galaxies of R = 27.0 ± 0.5, 
just 0.4 magnitudes fainter than the 3-cr limit of R = 26.6 
determined in §3.2. This suggests that these galaxies are a 
continuation of the (R — K) > 6.0 population that is de- 
tected in R, rather than a distinct population of galaxies 
with much more extreme colours. The three galaxies with 
the most extreme (R — K) colours are shown in Fig. 7. 

Pozzetti & Mannucci (2000) develop a method for 
breaking the "age-dust" degeneracy between passive and 
starburst EROs. They explore the (I — K)-(J — K) and 
(R — K)-(J — K) planes, proposing a classification scheme 
based on the steepness of the spectral break between optical 
and near-infrared bands at z — 1-2. We show the (R — K)- 
(J — K) plane of those EROs for which we have obtained 
J-band detections in Fig. 6b. We note that all of the mor- 
phologically classified "Compact" EROs (see §3.5) lie on the 
"Elliptical" side of the classification boundary. The situation 
is less clear for the "Irregular" galaxies, where the large pho- 
tometric uncertainties mean that we cannot reliably state 
that the majority of these galaxies lie on the "Starburst" 
side of the classification boundary. More accurate J-band 



photometry and spectroscopic identifications of the brighter 
members are required before firm conclusions may be drawn 
from this analysis. 

Whilst the foreground cluster lenses increase the 
chances of successful spectroscopic identifications, the same 
amplifying power allows us to observe yet fainter EROs 
(K ~ 20-21 in the image plane), for which spectroscopic 
observations will remain unfeasible. As a test of the ca- 
pabilities of photometric redshift measurements for EROs 
we use hyper-z (Bolzonella et al. 2000) to estimate pho- 
tometric redshifts for those EROs for which RIJK pho- 
tometry is available. We use template spectra correspond- 
ing to 51 different ages for a r = 0.1 Gyr p-model from 
Bruzual & Chariot (1993) to obtain a likelihood map in 
dust extinction-redshift (Ay-z) space for each galaxy. In 
most cases, valid solutions are possible for a wide range of 
redshifts, however the most likely solution is usually found 
around 0.8 < z < 1.5-2. 

Future spectroscopic observations (e.g. Smith et al. 
2001b) in conjunction with J- and //-band photometry of 
the whole sample will probe the redshift distribution and 
star formation histories of our ERO sample in significantly 
more detail. 

3.5 Morphological Diversity 

To gain further insight into the diversity of the ERO popu- 
lation we now attempt to quantify the morphologies of our 
sample. The optical imaging from HST has exquisite spa- 
tial resolution (~0.l"), however EROs are optically faint 
(R <; 24), and therefore suffer from low signal to noise in 
optical pass-bands. In contrast, EROs are typically well de- 
tected in the near-infrared and our A"-band imaging from 
UKIRT enjoys superb seeing (~0.5"). We therefore exploit 
both datasets in the following analysis. 

We first attempt to measure the scale-size of each mem- 
ber in our ERO sample. We fit a Moffatt profile to the 
/f-band image of each ERO from which we estimate the 
FWHM, and then crudely correct for the effects of seeing 
by subtracting in quadrature the FWHM of the PSF on the 
frame (Table 3). These measurements reveal no trend in in- 
trinsic FWHM versus colour or /f-band magnitude, most 
likely due to uncertainties in the individual measurements. 

We therefore attempt to measure the properties (scale- 
size and central concentration) of a typical ERO in ampli- 
fication corrected magnitude bins: K < 19, 19<ZT<20 
and K > 20, co-adding the Zf-band images of the EROs 
to obtain one representative galaxy per bin. Using r\ (0) = 
1(0) /(I)e where 1(0) is the surface brightness at radius 
and (I)e is the mean surface brightness within (Petrosian 
1976; Kron 1980), we measure the angular size So. 5 of the 
three co-added galaxies where rj (60.5) = 0.5 (Bershady et al. 
1998). We also measure the concentration of these composite 
galaxies, C v = F (< 0o.s) / F (< 1.5 So. 5), i.e. the ratio be- 
tween the flux within the radius So. 5 and that within 1.5 So. 5 
(Saracco et al. 1999). 

Despite the co-addition of the data in these broad mag- 
nitude bins, no discernible trend in scale-size or concentra- 
tion with /f-band magnitude or colour is found; for exam- 
ple, the concentration values differ by g 0.01 while the er- 
ror bars are typically Z 0.05. It therefore appears that, al- 
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Figure 5. The cumulative surface density of (a) (R — K) > 5.3 and (b) (R — K) > 6.0 EROs detected through our cluster lens sample 
after removal of gravitational amplification, as a function of limiting source-plane A"-band magnitude. The faintest two bins in both 
panels are plotted after allowing for the effects of incompleteness. We make no pass-band corrections between the different i\-band 
filters used in the three surveys as such corrections are ;£ 0.1 magnitudes (Cowie et al. 1994; Pcrsson ct al. 1998). The error bars are all 
Poissonian, and therefore generally underestimate the true error due to source clustering. We also plot the data of Daddi et al. (2000a) 
and Thompson et al. (1999), together with the PLE model predictions of Daddi et al. (2000b) and McCrackcn ct al. (2000) and the 
semi-analytic model prediction of Cole et al. (2000). Note that the McCrackcn et al. model predicts zero surface density of EROs with 
(R — K) > 6.0. The Daddi ct al. models assume 

(Q = 0.3, A = 0.7, H = SOkms^Mpc- 1 ), McCracken ct al. assume (f2o = 0.1, 
Ao = 0, Ho = 50kms~ 1 Mpc _1 ) and Cole et al. assume (f2o = 0.3, Aq = 0.7, Ho = 70kms~ 1 Mpc~ 1 ). Over the redshift range of interest 
(z ~l-2), the ACDM cosmology adopted by Daddi et al. and Cole et al. is very similar to McCracken et al.'s open cosmology (Peebles 
1984 & McCracken et al. 2000). These differences in cosmological model therefore do not affect the conclusions drawn in §4.1 and §4.2. 



though our K-b&nd data enjoys both high signal to noise 
and high resolution (0.5"), a detailed morphological analysis 
of these galaxies at near-infrared wavelengths is not possible. 
We therefore turn to the optical HST imaging, noting that 
at 2 ~ 1-2, the F702W filter samples the rest-frame ultra- 
violet. The following morphological classification is therefore 
sensitive to any ongoing unobscured star-formation in these 
galaxies, allowing us to identify easily starburst systems.. 

We classify our ERO sample on the basis of their ap- 
pearance in the HST F702W frames (Table 3): "C", com- 
pact galaxies; "I", irregular galaxies, e.g. disk-like, clumpy 
or interacting morphology; "F", faint galaxies which are 
either not or only just detected in the HST frames. The 
(R-K)> 5.3 sample comprises: 18% (C); 50% (I); 32% (F). 
Considering just the Compact and Irregular galaxies, we find 
that ~ 90% of the (R — K) > 6.0 population are irregular, 
compared with ~ 65% of those with 5.3 < (R — K) < 6.0. 
This is another hint that the (R — K) > 6.0 population may 
be dominated by distant dusty starbursts, although the large 
fraction of irregular 5.3 < (R — K) < 6.0 EROs, suggests 
that dusty starbursts may also make an important contri- 
bution to the less extreme ERO population. 

These results, specifically the large fraction of dusty 
starburst EROs in our sample, disagree with previous at- 
tempts to classify EROs on the basis of their morphology. 
For example, Moriondo et al. (2000) claim that only ~15- 
20% of EROs in their sample have an "irregular" morphol- 



ogy, while 50-80% are well fitted by an elliptical galaxy (r 1 / 4 
law) profile. As noted by these authors, their data is het- 
erogeneous, being drawn from the HST archive, and so the 
depth of their observations is not well defined, in contrast 
to our highly homogeneous data set. Another difference be- 
tween the two studies is that ~ 65% of Moriondo et al.'s 
EROs come from targeted searches in probable dense envi- 
ronments, such as high redshift clusters, and regions around 
radio galaxies and quasars, whereas we survey random fields 
in the z <; 1 universe. Moriondo et al.'s (2000) results may 
therefore be biased by their concentration on dense environ- 
ments. 

Finally, we illustrate the diversity of ERO morphol- 
ogy in Fig. 7, including candidate interacting galaxies and 
strongly lensed EROs. 

3.6 Clustering of EROs 

There is an order of magnitude variation in the number 
of EROs detected in each cluster field (Fig. 2 & Table 1), 
supporting previous claims that at least some component 
of the ERO population is strongly clustered (Daddi et al. 
2000a). However, to confidently identify two or more EROs 
as lying at the same redshift and therefore as being phys- 
ically associated with each other generally requires spectra 
of those galaxies. Here, we search our photometric catalogue 
for examples of two or more neighbouring EROs that ex- 
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Figure 6. a) (R — K)—(I — K) colour-colour diagram showing our sample of EROs with (R — K) > 5.3. The predicted evolution of an 
early type galaxy in the Zf = 2.5, r = 0.1 Gyr PLE model (Daddi et al. 2000b) is shown by the solid line. We mark redshift intervals of 
Az ~ 0.25 by stars, starting with z = 1 at the lower left. The excess of galaxies with colours redder than the elliptical galaxy evolutionary 
track suggests that the (R — K) > 6.0 ERO population may contain a significant fraction of dusty starburst galaxies. The size of the 
filled points scale with source plane if-band magnitude and open points represent galaxies for which only limits on both colours are 
available, b) (R — K)-( J — K) colour-colour diagram showing the EROs in our sample for which we have obtained J-band detections. 
The solid diagonal line shows (J — K) = 0.34(R — K) + 0.19 from Pozzctti & Mannucci (2000). The dashed lines mark our ERO selection 
criteria. The size of the points scale with source plane if -band magnitude and the symbols denote the morphological classification of 
each ERO from §3.5. 



hibit similar A'-band magnitudes and (R — K) colours. The 
premise being that such systems may represent high-redshift 
galaxy associations and that the similar colours and magni- 
tudes of such galaxies are analogous with the strong early- 
type galaxy sequences observed in lower redshift clusters and 
groups (e.g. Fig. 2&3a). We find two candidates for pairs 
of EROs, one each in A 963 and A 1835. 

EROJ101701+39034 & ERO <J '101703+3903 '4 - 
These two EROs (Table 3) have K ~ 19 and (R - K) ~ 6 
and lie north-west of the cD galaxy in A 963 in Fig. 1. They 
both have a regular 7f-band morphology and very faint, low- 
surface brightness R— and J-band morphologies. Within the 
photometric errors, the (R — K) and (J — K) colours of these 
two galaxies support the idea that these are two high redshift 
elliptical galaxies. However, ERO J101701+3903.4 is ~0.5- 
1.0 magnitude bluer than its neighbour in (I — K) (Fig. 5a 
& Table 3), indicating that it may be at a slightly lower 
redshift. 

ERO J 140057+0252.4 — This is one of two adjacent 
galaxies (separation ~3") in the field of A 1835. The other 
galaxy falls just below the (R — K) > 5.3 criterion (Fig. 2), 
and is therefore not included in the ERO sample. Both 
galaxies have featureless Tf-band morphologies, as does the 
ERO in the i?-band. The neighbour however has a very dif- 
fuse and low surface brightness i?-band morphology. 

We conclude that there are no unambiguous examples 
of ERO associations in our sample. 



3.7 Strongly Lensed EROs 

The amplifying power of the cluster lenses becomes very 
high within <; 30" of the centre of the cluster. This is the 
region of the image plane in which rare, highly magnified 
giant arcs are detected in some clusters (e.g. A 383, Smith 
et al. 2001a). In this same region of the image plane, we 
detect three strongly lensed EROs in our cluster sample, 
which we discuss below. 

EROJ003707+ 0909.4 & EROJ003707+0909.5 - 
Three images of a single background galaxy are detected 
in the core of A 68. The two brighter images are adjacent 
to the central galaxy and are shown in Fig. 7, while the 
third image (ERO J0037006+0909.1) lies ~ 20" south of the 
central galaxy and is considerably fainter. The optical mor- 
phology of these images reveals complex structure, contain- 
ing what appears to be five bright knots within each image. 
In contrast the 7\-band shows a bright centrally concen- 
trated source (K=17.2, 17.6 and 19.1 respectively). We ob- 
tain colours of (R-K)~6,(I-K)~5 and ( J - K) ~ 2.3 
for the central red region in both of the two brighter im- 
ages (Table 3). The morphology and position of these arcs 
relative to the central galaxy suggest that they are the im- 
ages of a background spiral galaxy, probably at a redshift of 
z ~ 0.5-1.0. 

We also perform photometry using concentric circular 
apertures of increasing diameter to estimate the colour gra- 
dient within the brightest image. The colour gradient indi- 
cates that this galaxy only has a colour of (R — K) > 5.3 
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Figure 7. Images of a selection of EROs from the sample, showing the if-band to the left and the F702W-band (smoothed with a 0.2" 
FWHM gaussian) to the right in each case. The top row shows the most extreme EROs with (R — K) 7, the next row (upper middle) 
shows examples of strongly lenscd EROs, the third row (lower middle) illustrates three of the more extended systems and the bottom 
row gives apparently interacting or morphologically complex systems. Each panel has North top and East left and is 10" square. 



for photometric apertures of g 3" in diameter. From our 
lens model we estimate that the amplification of the cen- 
tral region of the galaxy is ~ 15, implying that a source 
plane aperture of diameter <; 0.75" would be required to 
identify this as an ERO. This source would therefore have 
(R — K) < 5.3 (2" diameter aperture - see §3.1) and would 
not be detected as as ERO in the absence of the foreground 
cluster lens. 

EROJ133521+4100.4 — This ERO is ~ 20" north east 
of the central galaxy in A 1763 and is shown in Fig. 7. Its 
ellipticity (a / b ~ 5) and position angle with respect to the 
central galaxy indicate that it is probably strongly lensed. 
We do not however detect any counter images which im- 
plies that this galaxy is either singly imaged, or that the 
counter images fall below the detection threshold of our ob- 
servations. The arclet has K = 19.8 and (R — K) ~ 6.6 and 
we place lower limits on its (I — K) and (J — K) colours of 
3.5 and 1.5 respectively, suggesting that this may be a high 
redshift elliptical galaxy. 

ERO J024804- 0332.1 — This ERO was identified by 
Smith et al. (2001a) and lies ~ 25" south of the cD galaxy 
in A 383 and adjacent to a bright cluster elliptical. The el- 
lipticity (a / b ~7) and position angle (tangential to the cD) 
of this ERO supporting the interpretation of this image as a 



lensed background galaxy. We measure (R — K) ~ 6 and 
(I — K) £ 5, making this an extreme ERO under both 
(R — K) and (I — K) selection criteria. Smith et al. (2001a) 
estimated from their lens model that the object's redshift is 
U4. 



4 DISCUSSION AND CONCLUSIONS 

In this section we compare the observed number counts of 
EROs with the theoretical predictions of pure luminosity 
evolution (Daddi et al. 2000b; McCracken et al. 2000 - here- 
after MOO) and semi-analytic (Cole et al. 2000 - hereafter 
COO) models of galaxy formation, and then summarise our 
conclusions. Our primary aim when comparing our obser- 
vational results with the different model predictions is to 
investigate the suggestion that EROs comprise a mixture of 
evolved and dusty galaxy populations. 



4.1 Comparison with PLE Models 

In Fig. 5 we compare our observed number counts (§3.3) 
with Daddi et al.'s (2000b) PLE model down to K ~ 22. 
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This model attempts to describe EROs as a single popula- 
tion of evolved galaxies using two free parameters: Zf, the 
initial redshift of star formation in these galaxies and r, the 
c-folding time of the starburst. Daddi et al. adopt a Salpeter 
(1955) initial mass function (IMF), solar metallicity, no dust 
reddening and normalise the predicted number counts to the 
local luminosity function and rest frame colours of elliptical 
galaxies. Looking at the range of model parameters in this 
figure we see that at K ~ 21 the models diverge, enabling us 
to discriminate between them. For the (R — K) > 5.3 sam- 
ple, the models with r = 0.1 Gyr, Zf — 2.5 and r = 0.3 Gyr, 
Zf = 3.5, best match the data, with higher values of Zf pre- 
dicting too many EROs at K £ 20. The same models are 
also the best match with the (R — K) > 6.0 sample. 

However, as we show earlier, EROs appear not to be a 
single population of galaxies (§1), with the fraction of dusty 
starbursts possibly increasing towards redder optical/near- 
infrared colours and fainter if -band magnitudes (§3.3, §3.5). 
Agreement between the observations and these simple mod- 
els should therefore deteriorate, at the faint, red limit. Un- 
fortunately, comparison close to the faint limit is hampered 
by the cumulative nature of the data plotted in Fig. 5, fainter 
magnitudes being influenced by the brighter bins. Neverthe- 
less, Fig. 5b suggests that the models have more difficulty 
in correctly predicting the surface density of (R — K) > 6.0 
EROs, supporting the idea that this population contains a 
substantial fraction of dusty starbursts. We also note that 
although the best model lies within 1-a (Poissonian) of our 
data, the very small number of these redder EROs in our 
sample, particularly at K ~ 18-19, implies that the dis- 
agreement between Daddi et al.'s (2000a) data and this 
model at K ~ 18-20 is a better test of the model at 
these magnitudes. Daddi et al. (2000b) also noticed this dis- 
crepancy between their (R — K) > 6 data and the mod- 
els, and suggested that it was due to a possible deficit of 
(R - K) > 6.0 ellipticals at K g 19.5. 

In summary, a simple model that assumes EROs com- 
prise only passively evolving elliptical galaxies (with r = 
0.1 Gyr and Zf = 2.5) succeeds fairly well in predicting the 
number density of the (R—K) > 5.3 population. This model 
is also the best match to observations of the (R — K) > 6.0 
population, however the agreement is significantly worse, 
suggesting that the "single population" assumption may be 
a poor description of these more extreme EROs. 

We also compare our observations with a more general 
PLE model which attempts to describe the whole galaxy 
population, and not just passive ellipticals. This more de- 
tailed model (MOO) includes five galaxy populations (E/S0, 
Sab, Sbc, Scd, Sdm), each normalised to observed local 
galaxy parameters (i.e. luminosity function and rest frame 
colours). MOO also require their model to reproduce the 
shape and amplitude of Cowie et al.'s (1996) K -selected red- 
shift distribution, which contains very few K < 19 galaxies 
at z > 1. To achieve this, they adopt a dwarf dominated 
IMF (a; = 3) in order to reconcile their PLE prescription 
with Cowie et al.'s observational results. 

The MOO model (Fig. 5) under-predicts the surface den- 
sity of (R — K) > 5.3 EROs by approximately an order of 
magnitude and predicts none with (R — K) > 6.0. We sug- 
gest that this is probably caused by MOO requiring their 
model to fit Cowie et al.'s (1996) -ff-selected redshift distri- 
bution, which contains comparatively few galaxies at z > 1. 



However, Cowie et al. (1996) were concerned that their op- 
tical follow-up of a small-field Tf-selected sample might be 
incomplete for the reddest (i.e. optically faintest - R <; 24) 
galaxies, and hence for galaxies at z ^ 1. Recent wide field 
surveys have also discovered that EROs are strongly clus- 
tered (e.g. Daddi et al. 2000a: ~ 700arcmin 2 ). In contrast, 
Cowie et al.'s spectroscopic survey covered just 26.2 arcmin 2 , 
raising the possibility that their survey targeted an under- 
dense patch of sky. We therefore suggest that near-infrared 
spectroscopy of wide field 7f-selected samples is necessary 
before the redshift distribution of such galaxies beyond z ~ 1 
can be reliably quantified. 

4.2 Comparison with Semi-analytic Models 

We now return to the primary aim of our comparison with 
model predictions, that of investigating the suggestion that 
EROs comprise a mixture of evolved and dusty galaxy pop- 
ulations. 

We compare our ERO number counts with the predic- 
tions of the "reference model" from COO. This is a semi- 
analytic model which calculates the formation and evolution 
of galaxies in hierarchical clustering cosmologies, based on 
N-body simulations and simple parametrisation of physical 
processes. The model parameters are constrained by a num- 
ber of local galaxy properties including the ratio of elliptical 
to spiral galaxies, the metallicity of local L* ellipticals, B— 
and if -band luminosity functions, the fraction of gas in spi- 
ral and irregular galaxies and the size of galaxy disks. COO 
assume a Kennicutt (1993) IMF, however they find that the 
optical/near- infrared colours of galaxies are insensitive to 
this choice. They also include dust extinction using a Milky 
Way extinction curve, however as the authors point out, 
they do not allow for clumping of the dust and stars, nor 
do they include the effects of dust emission. We note that 
the model predicts a much weaker colour-magnitude corre- 
lation for cluster elliptical galaxies than that observed in 
the Coma cluster, the predicted colour of cluster ellipticals 
being ~0. 1-0.4 magnitudes bluer than the observations. 

The surface density of EROs predicted by the COO ref- 
erence model is shown in Fig. 5 (Baugh, private communi- 
cation). Despite its success in reproducing the properties of 
local galaxies, this model under-predicts the surface density 
of EROs with (R — K) > 5.3 and 6.0 by approximately an 
order of magnitude. It would obviously be inappropriate to 
modify this model to fit the ERO number counts if in do- 
ing so, it no longer agreed with the local observational con- 
straints mentioned above. Nevertheless, this disagreement 
with observations may point to important opportunities to 
improve our understanding of galaxy formation and hence 
improve COO's semi-analytic model. We therefore briefly con- 
sider where these opportunities may lie. 

As EROs appear to comprise galaxies containing both 
evolved and dust-reddened stellar populations (§1, §3.3 and 
§3.5), COO's under-prediction of EROs suggests that their 
reference model contains insufficient old stars and/or dust 
at z ~ 1-2. It therefore appears that their reference model 
does not produce enough stars and/or dust at high-redshift 
(z S> 2). One possible remedy would be to increase the 
fraction of stars formed in bursts, at the expense of quies- 
cent star formation. This should reduce the level of ongoing 
star formation at z > 1 (which would otherwise make these 
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galaxies too blue to be classed as EROs). Another possi- 
ble remedy would be to adopt a giant dominated IMF as 
this would increase the quantity of dust produced through 
the formation of a larger fraction of massive stars, thus both 
making starburst galaxies redder and allowing more systems 
to go through a starburst ERO phase without altering the 
properties of the local galaxy population with which the 
current model agrees. Introducing dust clumping may also 
provide a more realistic description of the spatial distribu- 
tion of dust in starburst systems. These options should be 
considered in future revisions of the semi-analytic models. 

Although the COO reference model fails to reproduce 
the observed number of EROs, this model does provide 
a physically motivated framework to qualitatively inter- 
pret the relationship between the various classes of ERO. 
In the model the redshifts of galaxies predicted to have 
5.3 < (R - K) < 6.0 and K £ 22 are z ~ 0.8-2.5, equally 
split between the z ~ 0.8-1.5 and z ~ 1.5-2.5 bins. The 
lower redshift bin is dominated by passively evolving galax- 
ies (~ 90% of the total) whose stars are sufficiently old to 
produce such red colours in the observed passbands, whilst 
the higher redshift bin is dominated by starburst galaxies 
(~ 95% of the total), whose (R — K) colours are reddened 
by dust. A similar analysis of the predicted (R — K) > 6.0 
population reveals that almost all of these more extreme 
EROs are predicted to lie at z ~ 1.5-2.5, with ~ 70% of 
these galaxies being dusty starbursts. The COO reference 
model therefore appears to qualitatively support the obser- 
vational evidence (§3) for both a predominance of distant 
dusty starburst galaxies in the (R — K) > 6.0 EROs, and 
passively evolving galaxies at z ~ 1 producing the bulk of 
the (R — K ) ~ 5.3-6.0 population. 

We now return to the break in the number counts 
of EROs identified in §3. Whilst the normalisation of the 
semi-analytic predictions falls an order of magnitude short 
of our observations, the qualitative properties of the pre- 
dicted ERO population offers a plausible explanation of 
the break in the count slopes. At magnitudes brighter than 
the break, the counts may be dominated by lower redshift 
(z ~ 0.8-1.5), predominantly passively evolving galaxies, 
whereas faint-ward of the break, the counts could be dom- 
inated by the distant (z ~ 1.5-2.5), dusty starburst galax- 
ies. We therefore speculate that the break may be due to a 
transition from an ERO population dominated by evolved 
galaxies at z ~ 1 (K <; 19.5) to one dominated by distant 
dusty starburst galaxies (that may have experienced a recent 
merger) at z ~ 2 (K £ 19.5). 

4.3 Conclusions 

We have undertaken a deep optical/near- infrared survey of 
10 massive cluster lenses at z ~ 0.2 using HST and UKIRT. 
We find 60 EROs with (R — K) > 5.3, of which 26 have 
(R — K) > 6.0 in a total image plane survey area of 49 
arcmin down to K = 20.6. 

We use detailed models of the cluster lenses (S01) to 
quantify the lens amplification and thus correct the ob- 
served number counts for the effects of gravitational lens- 
ing. After making these corrections, we estimate a surface 
density of 2.5 ± 0.4 (1.2 ± 0.3) arcmin -2 for EROs with 
(R - K) > 5.3 (6.0) at K < 21.6. Our results agree with 
previous shallower wide-field surveys at K g 19 and probe 



the number density of EROs with (R — K) > 5.3 and 6.0 
down to a source plane magnitude of K ~ 22. 

The number counts of both classes of ERO flatten sig- 
nificantly at magnitudes fainter than K ~ 19-20. We specu- 
late that this is due to a transition from an ERO population 
dominated by evolved galaxies at z ~ 1-2 (K g 19-20) to 
one dominated by distant dusty starburst galaxies at z > 1 
(K :> 19-20). Analysis of the (R-K)-(I-K) and (R-K)- 
(J — K) planes also suggests that the (R — K) > 6.0 popu- 
lation may contain a substantial fraction of dusty starburst 
galaxies. 

Approximately 50% of our sample contain morpholog- 
ical substructure including disk-like, clumpy or interacting 
morphologies. This is a larger fraction than found by previ- 
ous studies (e.g. Moriondo et al. 2000) which claimed that 
only ~15-20% of EROs have such morphologies. The dis- 
crepancy may be due to a bias towards dense environments 
and the heterogeneity of Moriondo et al.'s dataset, compared 
with our unbiased and more homogeneous dataset. 

We compare our observations with progressively more 
sophisticated models of galaxy formation, beginning with 
a two parameter (zf and r) PLE model that attempts to 
describe EROs as a single population of elliptical galax- 
ies (Daddi et al. 2000b). The model parameters which best 
match the observations are r = 0.1 Gyr and Zf = 2.5, ruling 
out the very high formation redshifts (z/ ~ 10) that were al- 
lowed by Daddi et al.'s (2000a) shallower observations. How- 
ever, this single population model matches the (R—K) > 6.0 
EROs significantly worse than the (R—K) > 5.3 EROs, sup- 
porting the idea that the more extreme population contains 
a large fraction of distant dusty starbursts in addition to the 
elliptical galaxies contained within this model. 

We then compared our observations with PLE mod- 
els that attempt to describe the whole galaxy population, 
and not just passive ellipticals (MOO). These models under- 
predict the surface density of (R — K) > 5.3 EROs by ap- 
proximately an order of magnitude and predict none with 
(R — K) > 6.0. This deficit of EROs is probably caused 
by MOO requiring their model to fit the median redshift of 
Cowie et al.'s (1996) _R"-selected redshift distribution. We 
suggest that this confirms Cowie et al.'s (1996) concern that 
their optical follow-up of a small-field if-selected sample 
might be incomplete for the reddest (i.e. optically faintest - 
R Z 24) galaxies, and hence for galaxies at z £ 1. It there- 
fore appears that near-infrared spectroscopy of wide field 
if-selected samples is necessary before the redshift distri- 
bution of galaxies at z Z 1 can be reliably quantified. 

Finally, we compare the observed number density of 
EROs with the semi-analytic predictions from the reference 
model of COO. This semi-analytic model under-predicts the 
number density of EROs at K ~ 18-22 by an order of magni- 
tude, indicating that the current generation of semi-analytic 
models may not produce sufficient stars and/or dust at high 
redshift (z > 2). However, as the COO reference model is 
physically well motivated, we look at the predicted proper- 
ties of the ERO population, finding that their redshift distri- 
bution and the split between passive and dusty EROs appear 
to support our interpretation of the break in the slope of the 
number counts. 

This is the first survey to exploit massive foreground 
galaxy clusters to amplify the flux of background EROs. 
The sample constructed from our deep, high resolution K- 
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and i?-band observations is therefore ideally suited to near- 
infrared spectroscopic follow-up on 10-m class telescopes 
(e.g. Smith et al. 2001b), the lens amplification allowing us 
to probe ~ 1 magnitude beyond the normal reach of such ob- 
servations. The primary goal of our future programme will 
be to probe the properties of EROs beyond the break in 
their number counts at K ~ 19-20 and to identify the cause 
of this feature. 
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Table 3. Catalogue of EROs 



Source 



a,8 
(J2000) 



K 



(R-K) 



(I -K) (J - K) 



FWHM 

(") 



Morphology 

(a) 



A 68 






















ERO J003703+0909.7 


00 37 03.30 


+09 09 44.4 


20.10 + 0.16 


> 6.39 


> 4.40 


> 


1.42 


1.23 


F 




ERO J003706+0909.1 


00 3706.10 


+09 09 08.7 


19.10 + 0.01 


> 5.70 


> 3.86 


> 


2.06 


0.92 


I, 


§3.6 


ERO J003707+0909.4 


00 3707.25 


+09 09 23.8 


17.62 + 0.03 


6.21 + 0.22 


5.25 ± 0.35 


2.32 


±0.18 


0.62 


I, 


§3.6 


ERO J003707+0909.5 
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+09 09 28.4 


17.16 + 0.03 


6.46 + 0.26 


5.41 ± 0.30 


2.29 


±0.15 


0.66 


I, 


§3.6 


ERO J003709+0909.6 


00 3708.84 


+09 09 37.0 


19.45 + 0.10 


6.19 + 0.20 


> 5.33 


> 


2.07 


0.25 


F 




ERO J003709+0911.0 


00 37 09.03 
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4.00 ± 0.03 


2.00 
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I, 


Clumpy 
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r\r\ TV C\C\ c T 
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0.18 
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-13 36 28.4 


18.96 + 0.09 


6.19 + 0.13 


4.74 + 0.10 






0.36 


F 




A 267 






















ERO J015238+0101.0 


0152 38.43 


+0102 02.8 


17.76 + 0.02 


5.51 + 0.05 


4.06 ± 0.05 


2.10 


±0.20 


0.42 


I, 


Disk? 


ERO J015238+0101.9 


0152 37.60 


+010151.6 


18.73 + 0.05 


6.15 + 0.12 


> 5.48 


> 


1.85 


0.67 


C 




ERO J015240+0101.2 


0152 40.08 


+010108.8 


18.74 + 0.05 


5.67 + 0.08 


4.52 ± 0.15 


1.96 


±0.30 


0.35 


C 




ERO J015240+0101.6 


0152 39.55 


+0100 38.4 


20.37 + 0.17 


> 6.65 


> 4.03 


> 


0.20 


0.57 


F 




ERO J015240+0101.7 


0152 39.88 


+010139.3 


20.32 + 0.15 


5.93 + 0.35 


4.09 ± 0.33 


> 


0.26 


0.34 


C 




ERO J015240+0101.8 


0152 39.82 


+010150.0 


19.51 + 0.10 


6.09 + 0.21 


> 4.59 


> 


1.07 


0.53 


F 




ERO J015241+0101.1 


015241.25 


+0102 07.4 


19.14 + 0.09 


6.22 + 0.21 


> 4.85 


> 


1.44 


0.68 


F 




ERO J015243+0101.5 


0152 42.80 


+010127.0 


19.93 + 0.11 


5.51 + 0.34 


4.78 ± 0.33 


> 


0.65 


1.30 


C 




ERO J015245+0101.8 


0152 44.53 


+0100 47.9 


20.04 + 0.16 


5.32 + 0.20 


> 4.21 


> 


0.54 


0.60 


F 




ERO J015247+0101.1 


0152 47.28 


+010104.8 


18.48 + 0.03 


5.72 + 0.07 


4.39 ± 0.07 


1.86 


±0.25 


0.29 


I 




ERO J015249+0100.6 


0152 48.93 


+0100 38.1 


20.58 + 0.16 


> 5.89 


> 4.65 


> 


0.03 


0.51 


F 




ERO J015249+0101.6 


015249.10 


+0100 36.8 


20.05 + 0.11 


5.51 + 0.24 


> 4.14 


> 


0.53 


0.52 


F 




A 383 






















ERO J024804— 0332.1 


02 48 04.12 


-03 32 06.4 


19.73 + 0.07 


6.22 + 0.28 


> 4.53 






0.29 


C 


§3.6 


ERO J024805-0330.0 


02 48 04.55 


-03 29 57.3 


19.59 + 0.09 


6.27 + 0.17 


4.51 ± 0.25 






0.46 


I 




ERO J024805— 0330.2 


02 48 04.56 


-03 3011.0 


18.12 + 0.02 


5.41 + 0.04 


3.29 ± 0.03 






0.00 


C 


Star? 


ERO J024805-0330.3 


02 48 04.82 


-03 3018.0 


18.06 + 0.04 


5.61 + 0.06 


4.39 ± 0.07 






0.47 


I 




ERO J024805-0330.3 


02 48 04.90 


-03 30 22.7 


20.06 + 0.11 


6.31 + 0.33 


> 4.16 






0.36 


c 




ERO J024805— 0330.4 


02 48 04.89 


-03 30 24.7 


19.18 + 0.08 


6.08 + 0.15 


4.64 ± 0.25 






0.33 


I 




ERO J024806-0331.5 


02 48 06.44 


-03 3127.3 


18.27 + 0.05 


6.63 ± 0.30 


4.74 + 0.15 






0.54 


I 




ERO J024808-0331.8 


02 48 07.99 


-03 3147.5 


19.44 + 0.10 


6.55 + 0.28 


> 4.50 






0.63 


I 




ERO J024809— 0332.9 


02 48 08.79 


-03 32 56.4 


19.47 + 0.11 


5.56 + 0.17 


4.80 + 0.33 






0.62 


F 




ERO J024809-0333.0 


02 48 08.58 


-03 33 00.0 


19.00 + 0.05 


5.77 + 0.17 


4.97 + 0.25 






0.61 


I, 


Merging? 


ERO J024810— 0332.4 


02 4810.02 


-03 32 25.8 


17.57 + 0.02 


5.39 ± 0.04 


4.08 ± 0.05 






0.61 


I 





A 773 
ERO J091802+5143.1 

A 963 



0918 02.50 +5143 06.4 20.23 + 0.11 5.59 + 0.18 



> 3.73 



> 1.23 



0.34 



ERO J101706- 
ERO J101707- 
ERO J101707- 



3901 



I, Disk 



4 


101701.42 


+39 03 24.6 


19 


32 + 


05 


6 


09 ± 0.25 


3.17 ±0.15 


> 


2.13 





47 


1 


4 


1017 02.52 


+39 03 22.0 


18 


84 + 


04 


6 


03 + 0.17 


4.13 ± 0.25 


2.11 


±0.11 





58 


F 


7 


1017 05.69 


+39 0141.6 


20 


53 + 


15 


6 


38 ±0.34 


> 2.94 


> 


0.80 





00 


F 





1017 05.79 


+39 02 02.3 


20 


16 + 


10 


5 


55 ±0.19 


> 3.13 


> 


1.04 





57 


F 


4 


101706.61 


+39 0122.1 


19 


51 + 


08 


5 


78 ±0.13 


> 3.86 


> 


1.71 





11 


I 


1 


101706.58 


+39 02 08.3 


19 


74 + 


09 




> 6.90 


> 3.49 


> 


1.27 





21 


I, 


4 


1017 07.09 


+39 02 25.0 


20 


37 + 


09 


5 


37 ±0.17 


> 3.06 


> 


0.96 





28 


I 



Clumpy 



(a) "C", Compact; "I", Irregular; "F", Faint morphological classification, as discussed in §3. 
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Table 3. (continued) Catalogue of EROs 



Source 



a, S 
(J2000) 



K 



(R -K) (I- K) 



(J - K) FWHM Morphology 

(") (a) 



A 963 (continued) 

ERO J101707+3902.9 
ERO J101709+3903.5 
ERO J101710+3902.2 
ERO J101710+3903.0 
ERO J101710+3903.4 

A 1763 
ERO J133511+4100.3 
ERO J133517+4058.7 
ERO J133518+4058.8 
ERO J133518+4059.4 
ERO J133521+4100.4 

A 1835 
ERO J140057+0252.4 

A 2219 
ERO J164023+4644.0 
ERO J164024+4643.7 



101707.43 +39 02 52.2 
1017 09.30 +39 03 31.4 
1017 09.66 +39 0210.0 
1017 09.61 +39 03 03.6 
101709.94 +39 03 23.3 



13 3511.46 +410015.5 
13 3516.77 +40 58 43.0 
13 3517.82 +40 58 45.9 
13 3518.30 +40 59 25.8 
13 35 21.16 +4100 24.4 



19.14 + 0.04 

17.15 + 0.01 
19.60 + 0.08 
19.49 + 0.07 
19.66 + 0.09 

18.92 + 0.06 
20.11 + 0.15 
20.47 + 0.16 
17.65 + 0.02 
19.84 + 0.10 



5.63 + 0.14 
5.61 + 0.04 

> 7.28 
5.99 + 0.17 
6.46 + 0.21 

5.57 + 0.09 
5.44 + 0.13 
5.42 + 0.20 
6.08 + 0.06 
6.57 + 0.51 



> 4.29 
3.97 + 0.05 

> 3.90 

> 3.70 

> 3.60 

> 3.81 

> 3.26 

> 2.70 
4.83 ± 0.25 

> 3.05 



1.94 + 0.10 
1.94 + 0.05 

> 1.65 

> 1.55 

> 1.47 



> 1.57 

> 0.98 
2.32 + 0.20 

> 1.05 



14 00 57.13 +02 52 26.4 18.01 + 0.02 5.33 + 0.08 3.92 + 0.15 1.85 + 0.07 



16 40 23.05 +46 44 02.3 
16 40 23.95 +46 43 42.1 



17.51+0.01 
19.64 ± 0.05 



5.90 + 0.04 
5.71 + 0.12 



4.51 ± 0.08 
4.07 ± 0.30 



2.28 + 0.10 
> 1.76 



0.46 
0.56 
0.64 
1.20 
0.85 

1.25 
0.63 
0.64 
0.76 
1.24 

0.53 

0.51 
0.68 



I, Disk? 

C 

I 

I 

F 

I 
I 
I 
I 

F, §3.6 
C 

I, Disk 
I 



(a) "C", Compact; "I", Irregular; "F", Faint morphological classification, as discussed in §3. 
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